Critical behavior of the Random-Field Ising model at and beyond the Upper Critical 

Dimension 
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The disorder-driven phase transition of the RFIM is observed using exact ground-state computer 
simulations for hyper cubic lattices in d = 5, 6, 7 dimensions. Finite-size scaling analyses are used to 
calculate the critical point and the critical exponents of the specific heat, magnetization, susceptibil- 
ity and of the correlation length. For dimensions d — 6,7 which are larger or equal to the assumed up- 
per critical dimension, d u = 6, mean-field behaviour is found, i.e. a = 0, (3 — 1/2, 7=1, v = 1/2. 
For the analysis of the numerical data, it appears to be necessary to include recently proposed 
corrections to scaling at and beyond the upper critical dimension. 

PACS numbers: 64.60. De, 75.10.Nr, 75.40.-s,75.50.Lk 

Keywords: random-field Ising model, upper critical dimension, correction to scaling, mean field 



I. INTRODUCTION 

The random-field Ising model (RFIM)i is a prototypi- 
cal model for magnetic systems exhibiting quenched dis- 
order. For d < 3 and higher dimensions^, the model ex- 
hibits a phase transition: For low temperatures and/or 
weak disorder the coupling dominates and the system is 
long range ordered. Increasing the temperature and/or 
the disorder to a critical value, the RFIM with a Gaus- 
sian distribution of the disorder undergoes a second order 
phase transition to a paramagnetic phase. At dimensions 
below the lower critical dimension d < di — 2 the critical 
disorder strength is zero and no ferromagnetic long-range 
order— ^ occurs for finite values of the randomness. In 
d = 3 the RFIM is believed^ to describe the behavior 
of experimentally accessible antiferromagnets well. To 
investigate d = 3, 4 dimensional RFIM systems, many 
studies were published, see, e.g., Refs£~— and critical 
exponents for the magnetization (3 and the susceptibility 
7 are known with good accuracy. Results for the specific 
heat exponent a are diverse, varying from quite negative 
to small positive values. As a result some values support 
the scaling relation a + 2(3 + 7 = 2. Also a different 
scaling relation a + 2(3 + 7 = 1 has been proposed 7 . 

The RFIM has been studied by mean-field techniques 
as welU^r— . The upper critical dimension d u , above 
which the RFIM shows dimension-independent mean- 
field behavior, has been predicted to be d u — 6 1 . The 
critical exponents are believed to hold the mean field 
values of the pure Ising model. To the knowledge of the 
authors, no numerical simulations have been performed 
to confirm the value of the upper critical dimension of 
the RFIM or to study its critical behavior close to or at 
d u . Also corrections to scaling are still unknown in this 
region. In this paper we present our numerical results of 
the critical exponents around the upper critical dimen- 
sion in d — 5, 6, 7. Therefore exact ground states (GS) 
of large instances are calculated using a mapping to a 
combinatorial optimization problem. The main result is 
to confirm d u = 6 and all exponents for d u > 6 to be 
indeed compatible with the mean-field values of the pure 



Ising ferromagnet, if corrections to scaling are taken into 
account. A summary of this work can be found at the 
papercore database^. 

The RFIM consists of a hyper cubic lattice of Ising 
spins Si = ±1. The Hamiltonian of the RFIM is given 
by 



u 



(i) 
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J denotes the ferromagnetic coupling constant between 
two adjacent spins. denotes to sum over next 

neighboured spins only. Each spin is exposed to a net 
random-field hi + H which consists of two contributions: 
hi = hsi is the quenched local random-field with a dis- 
order strength h. The £j are distributed according to a 
Gaussian with zero mean and unity width. H indicates 
the strength of a homogeneous magnetic field, which is 
used to determine the susceptibility, but is set to H = 
elsewhere. 

The paper is organized as follows: In the next section 
we briefly outline the numerical approach we have used. 
Also we state in detail all quantities we measured and 
the finite-size scaling (FSS) approaches we have used to 
analyze the data. In section three we present our main 
results. In the last section, we discuss our work and con- 
clude with a summary. 



II. METHOD AND MEASURED QUANTITIES 

From mean-field theory, the phase transition of the 
Gaussian RFIM is known to be second order along the 
whole phase boundary-^. Furthermore, renormalisation 
group (RNG) theory brought up the existence of three 
fixed points for the renormalization group flow in the 
temperature/disorder phase space of the RFIM 20 . There 
exists a stable fixed point at T = h = 0, an unsta- 
ble fixed point at T = T c , h = and a saddle-point 
at T = 0, h = h c , see Fig. Q] The saddle-point is sta- 
ble against changes of the temperature and unstable to 
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pertubations of the disorder h. A crossing of the phase 
boundary at finite temperatures (T > 0) leads to the 
same critical behaviour as a crossing at T = by varying 
h does. This allows us to focus on ground states only, 




FIG. 1: Schematic phase diagram and renormalization group 
flow of the RFIM. 



which is very convenient, especially from the computa- 
tional point of view. We calculate the GS by mapping 
each realization ({hi},H) of the net random-field to a 
graph with suitable edge capacities. For this graph we 
apply a sophisticated maximum flow/minimum cut algo- 
rithm known from algorithmic graph theory to calculate 
the minimum cut^ir— . 

The cut separates the graph into two parts. They di- 
rectly correspond to the GS spin configuration {Si} of 
that specific realization of the net disorder. This ap- 
proach allows us to obtain exact ground states of large 
systems with a linear lattice lengths up to I = 14 in 
rf = 5, L = 9 in d = 6 and L — 8 in d — 7. In contrast, 
Monte Carlo would need an unresonable amount of CPU 
time for the same data quality, due to long equilibration 
times. 

From a GS spin configuration, we calculated directly 
the quantities of interest, which are magnetization per 
spin (JV is the number of spins) 
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and the bond energy per spin 
Ej 



N 51 SiS 3 



(2) 



(3) 



From these individual values, we obtained averaged quan- 
tities like the average magnetization m = [M]h, where 
[•••]/, denotes the disorder average for a given value of h. 
We also calculated the zero-temperature susceptibility 



X{h) = 



dM({h t },H) 



dH 



H=0 



(4) 



The actual value of the susceptibility is obtained as lin- 
ear response of the single-disorder realization magneti- 
zation to small homogeneous magnetic fields H. There- 
fore, we apply small homogeneous fields at equidistant 
values H\ , 2H\ , 3Hi to each realization of the disorder 
and (automatically) fit parabolas as function of H to the 
obtained single realization magnetizations M(H = 0, h), 
M(H = H u h), M (H = 2H u h) and M(H = 3Hi,h). 
The average of all linear coefficients of the parabolas at 
a particular h then is the susceptibility x(h) ( see Figure 
[2]). H\ is chosen depending on the dimension and the 
size of the system, see Sec. Mil 

We also calculate a specific heat like quantity C(h) as 
numerical derivative of the bond energy per spin Ej (see^ 
for details) which is defined as 



C(h) 



[9Ej(h)] t 
dh 



(5) 



From the magnetization at H = we deduce the Binder 
cumulant^, 



g(h,L) 



[M 4 



[M 2 
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We are interested in the infinite-size limit and the crit- 
ical behaviour of these quantities, especially at the crit- 
ical point. Our analysis proceeded in the following way: 
For each quantity we start at d = 5, where standard 
FSS arguments according to2£ are valid, continuing with 
d — 7, in which a correction-to-scalilng theory exists. Af- 
ter that we take care of its scaling in d = 6. We started 
at the "easy" quantities, as the Binder parameter and 
proceeded with the magnetization, susceptibility and at 
last the specific heat. 

In d = 5, standard FSS says, that finite-size effects 
enter via the lattice length to correlation length ratio 
i/Coo- Since the correlation length at the critical point 
diverges as £oo ~ \h — h c \~ u , one gets £/£oo ~ L\h — 
h c \ u . The same is valid for the Binder cumulant, g(h, L). 
Taking the argument to the power of one expects 
the Binder cumulant to scale as 

g(h,L) = g([h-h c )L l ' v ) (d<d u ), (7) 

= g(xW(L,h)) 

Equality is a result of the dimensionless definition of the 
Binder cumulant. Using a properly chosen /i-axis rescal- 
ing function X( d \L, h), the Binder cumulants for all sys- 
tem sizes belonging to the same dimension collapse on 
one single (unknown) master curve. Hence, finding the 
critical exponent v is now left as a problem with one de- 
gree of freedom: it is the value of v where the collapse of 
the rescaled data points from different system sizes onto 
one curve is best. 

The argument of the rescaling function in Eq. ([7]) van- 
ishes at the critical point h = h c . Thus the curves of 
the Binder cumulant do intersect. But for d > d u strong 
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corrections to scaling arise and the point of intersection 
for different system sizes drifts systematically. This phe- 
nomenum was discussed by Jones and Young2I for the 
pure Ising model. They found that for d > d u the lattice 
length has to be replaced by a larger one, i.e., 



a%L ln 1 ^" L d = 6 
L d/d u d > d 



(8) 



Thus, the scaling exponents of the correlation length de- 
pend on the dimension of the considered setup. In d = 7 
~ d u /d = 0.5-6/7 « 0.4286. For d = d u , 



it is v 



,MF 



the standard finite-size scaling arguments have to be 
adapted. Finding a correct form for X^ du '(L,h) was a 
tedious task. The first approaches to find the scaling 
form at the upper critical dimension began in the early 
1980's, when Brezin studied the </> 4 -model (d u — 4) on 
hypercubic d-dimensional lattices^. Using RNG calcu- 
lations he found multiplicative logarithmic corrections to 
scaling at the upper critical dimension. 

Much later, in 1997 Luijten and Blbtc 29 published a 
scaling form for 0(n) spin models with tunable ferromag- 
netic long-range interactions and tunable upper critical 
dimension. Their RNG calculations differ from Brezin's 
result by 2nd order logarithmic corrections. For the Ising 
model with long-range interactions at d = d u = 2 their 
scaling form was numerically verified by Griineberg and 
Hucht 30 .. But since the model contains no disorder we 
adapted some parameters and use the ansatz 

X^\h.L) = (/i-/i c )L 1 / ,/MF log(L) 1/6 + 61og(L) e (9) 

for FSS at d = d u = 6. 

But now, let us return to the other quantities of in- 
terest and lets discuss the scaling assumptions for them. 
The next less easy quantity is the magnetization. From 
scaling theory we expect the magnetization in d — 5 to 
scale as 

m{h,L) = L-^' v m(x^\L,h)\ (10) 
= L-P/ u m Uh-h c )L 1 ' v ^ . (11) 

so, replacing X^ H> X^ and L M> L should result in 
a data collapse for the higher dimensional setups. In 
the same way, a collapse of the susceptibility should be 
achieved, using 

X (h,L)=fr/»x(x W (L,hj) (12) 

For the scaling behavior of the specific heat like quantity 
we use 



C(h, L) = L- a '»C (X^(L, hj) , 



(13) 



since we are not aware of a suitable, analytically sound 
scaling form. 

We observe the behaviour of those quantities as func- 
tions for different linear lattice lengths L to apply finite- 
size scaling techniques. While the error-estimates of the 
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FIG. 2: (color online) Approximation of the susceptibility in 
d = 6 by fitting a parabola to different homogeneous magnetic 
fields using H\ — 0.036. For this realization we obtain via the 
coefficient of the linear term Xsingie ~ 5.2. 



observables are calculated from 200 bootstrap samples 3 -! 
the finite-size scaling analysis is done using two different 
tools from 3 - 2 - and 3 - 3 -. The former is a command line tool 
based on a simplex algorithm written in python, while 
the latter is a hand operated rescaling plot program with 
direct visualisation written in C. The simplex algorithm 
needs about 30 data points near the critical point to give 
reasonable results. As our data points for larger system 
sizes are relatively sparse compared to 30 but very ac- 
curate, we created artificial data points in that region. 
We did this by fitting cubic splines to our data points 
near the critical point. The splines are used to calculate 
a dense set of artificial data points to reach the necessary 
number to serve the simplex algorithm to create a reason- 
able data collapse. These data points are just auxiliary 
points, which means that we do not show any artificial 
data points in our final data-collapse plots. 



III. RESULTS 

We have performed exact GS calculations for the Gaus- 
sian RFIM in dimensions d — 5, 6 and 7. To perform a 
finite-size scaling analyses, we considered different sys- 
tem sizes up to L = 14 (d = 5), L = 9 (d — 6) and L = 8 
(d= 7). This means, we obtained exact GSs for systems 
with up to 2 x 10 6 spins. Each time, we averaged all re- 
sults over many realizations and studied each realization 
for several values of the disorder strength h. Note that for 
the calculation of the susceptibility, each time four GSs 
subject to the homogeneous fields H = Q,Hi,2Hi,3Hi 
were necessary. The homogeneous fields are chosen such 
that m(H) stays in the parabolic domain, see also the 
discussion in 9 -. For details see Tab. Q] Note that the 
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d = 


5 




d = 


6 




d = 


7 


L 


^real 


H l 


L 




#i 


L 




#i 








3 


60000 


0.0120 


3 


100000 


0.00350 


4 


80000 


0.00200 


4 


80000 


0.0050 


4 


80000 


0.00150 


5 


160000 


0.00075 


5 


30000 


0.0025 


5 


5500 


0.00075 


6 


39500 


0.00075 


6 


15000 


0.0015 


6 


2400 


0.00035 


7 


72000 


0.00050 


7 


7000 


0.0010 


7 


8000 


0.00025 


8 


40000 


0.00030 


8 


3500 


0.0006 


8 


5000 


0.00010 


9 


10000 


0.00100 


9 


1000 


0.0004 








10 


4400 


0.00020 














11 


2000 


0.00010 














12 


6800 


0.00030 














13 


3200 


0.00030 














14 


5000 


0.00025 















TABLE I: For each dimension d and linear size L, the max- 
imum number n rea i of realizations and the value Hi of the 
smallest homogeneous field used to obtain the susceptibility 
are listed. 



values of H± for d — 5 were chosen accidentally non- 
monotonously, but they are anyway small enough such 
that m(H) is parabolic, regardless of the system size. 

To get a first impression of the nature and location 
of the phase transition, we start by showing results for 
the Binder cumulant Eq. (J6j) . As mentioned above, for 
second-order phase transitions, the Binder cumulants ob- 
tained for different system sizes at or close to a critical 
point will intersect. The reader may get a rough guess of 
the critical disorder strength from the point of intersec- 
tion observing Fig. [3] Based on the data of the largest 
system sizes, our first estimates are — 6.02(1) in 

d = 5, hi e) = 7.78(1) in d = 6 and h ( J } = 9.48(5) in 
d = 7. The intersection is rather clear for d = 5. Using 
an appropriate magnification, it can be seen in the inset 
of Fig. [SJ that for d = 6 the points of intersection vary 
systematically. A similar also systematic drift may be 
observed for d — 7, but due to larger error bars, we do 
not quantify them. A more reliable way to determine the 
critical point is to apply finite-size scaling theory and col- 
lapse the data to a master curve. For d = 5 we follow the 
scaling ansatz according to Eq. ([7]). Our data collapse 



best is for hf ] = 6.0157(10) and v 



(5) 



d = 7 we used v — > v-j± obtaining h, 



0.626(10). For 
9.4889(1) and 

j/( 7 ' = 0.49(2) which is compatible with the mean-field 
result i^mf = 0.5. The collapsed data curves are shown in 
Fig. HI The data collapses look quite good for larger sys- 
tem sizes at d = 5, 7. When including the smallest sizes, 
the quality of the data collapses decreases. For d = 7 
it improved a bit by including some higher order correc- 
tions obtained by trial and error. Since we have no hint 
on the precise form of the corrections, we do not go into 
details and use the results from the high-quality collapse 
when omitting L = 3. 

In contrast to this, the data collapse for the Binder 
parameter for d = 6 looks only good for all system sizes 
when applying a scaling ansatz which includes logarith- 
mic corrections according to Eq. ^ . The best collapse is 




FIG. 3: (color online) Binder cumulant in d — 5, 6, 7 dimen- 
sions for different lattice sizes L = 4 ... 14, L = 3 ... 9 and 
L — 3 ... 8, respectively . The lines are guides to the eyes 
only. 
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FIG. 4: (color online) Collapsed Binder cumulants in d = 7 
and d = 5 (inset) for several system sizes. 



achieved using h { c 6) = 7.7859(1), = 0.51(5), c = 1/6, 
b = 3.403, and e = -3.0 (see Fig.[5j). 

The rescaling of the disorder axis according to 
X^ d \h,L) should be valid for all other quantities, so we 
will keep it fixed from here on. This allows us to de- 
termine the magnetization exponent /3 via rescaling the 
magnetization values. Using the same methods as above, 
we obtain in five dimensions (3^ = 0.255(10) for all sys- 
tem sizes. The collapsed data curves are shown in Fig. 

El 

More interesting are the cases for d = 6,7. Eq. 
© tells us to scale the magnetization in d = 7 with 
f3 /v d/d u = 1.16, which gives the best collapse for 
our data, but for the four largest system sizes only (see 
inset Fig. [7]). Due to the large deviation of the magnetiza- 
tion of the smaller system sizes, we tried to include some 



5 




FIG. 5: (color online) Collapsed Binder cumulants for d = 6 
(main plot). The inset shows the critical region where the es- 
timated curves of the Binder parameter intersect. The dotted 
black line corresponds to the value h = 7.7859. 




FIG. 6: (color online) Magnetization (main plot) and col- 
lapsed magnetisation (inset) in d = 5 for several system sizes. 



higher order terms. We achieved some collapses includ- 
ing these curves, but a definite form of the corrections is 
unclear. 

We also achieve good results in d — 6, see Fig. [JJ using 
a scaling form as 



rh(h) -m{h)L p '"h) h {L), 



(14) 



where = 0.50(1) and h = -0.33(1). It is also 
possible to collapse the data using a pure power law, 
m(h) ~ m(h)L om( -^ =>■ f3* = 0.40(1) utilizing the same 
disorder rescaling to achieve the same collapse quality. 
Nevertheless, the resulting value f3* is not compatible 
with the mean-field value /3mf = 0.5, see Tab. IIVI This 
indicates that indeed also here logarithmic corrections 
are important. 

We now turn to the susceptibility, as defined in Eq. (0} . 
It has a very clear peak at the critical point, which can 
be seen for example in Fig. |SJ Taking the peaks as a char- 



FIG. 7: (color online) Collapsed magnetisation in d — 6 for 
several system sizes and d = 7 (inset). 



acteristic property of the susceptibility we approximate 
them with parabolas. 

We use the maxima of the obtained parabolas for a 
FSS analysis. Following Jones and YoungSi we use the 
pure lattice length L in d = 5 and a stretched lattice 
length L for the larger dimensions, i.e. logarithmically 
enlarged L 6 = a 6 L ln 1 ^" L for d — 6 and L 7 = a 7 L d / du 
for d = 7 and try a pure power law ansatz according to 



Xmax(i) = s L l/u 



(15) 




FIG. 8: (color online) Susceptibility in d ■ 
the parabolas fitted to the peaks. 
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7. The inset shows 



For d — 5 we leave out the smallest systems for our 
analysis, since they do need unknown small-system size 
correction terms. The results of the fits are given in Tab. 
HB The original data points together with the fits are 
shown in Fig. |H1 

Since the values obtained for 7 are slightly larger than 
the mean-field values we check the results collapsing the 



6 



d 


s •y/v 


7 


5 


0.036(9) 2.27(11) 


1.42(7) 


6 


0.064(2) 2.13(2) 


1.07(3) 


7 


0.050(2) 2.24(2) 


1.06(7) 



TABLE II: Resulting fit parameters from fitting the peak 
height Xmax of the susceptibilities as a function of system 
size L according to Eq. ([15)1 . 




FIG. 9: (color online) Log-log plot of the maxima of the sus- 
ceptibilities as function of the lattice length for different di- 
mensions and their fitted power law model. The lines show 
the fits according to Eq. (|15|) . 



data. The finite-size behaviour of the disorder is taken 
from the results of the Binder-parameter analysis and 
stays unchanged from here on. The collapse of the sus- 
ceptibility of the RFIM at zero temperature at d = 6 
is shown in Fig. [10] using the fixed value v — 0.51(5), 

7 = 1.07(5) and — 7.7859. Given the error bars 
of the data points, the data collapse appears to be fair. 
The next quantity we present here is the specific heat 
as defined in Eq. ([5]). The curves of the specific heat 
for all system sizes in d = 5 are shown in Fig. 111! One 
can clearly observe a maximum, which lies slightly be- 
low h c . We determined for all dimensions and system 
sizes the positions and heights of the peaks by fitting 
parabolas to the peak locations. For the very large sys- 
tem, as d = 6, L = 9 and d — 7, L = 8 where the peaks 
of the specific heat are too ragged for a decent fit, we 
draw an ellipsoid around its most likely position, taking 
its radii as error bars. From the peak heights we observe 
the finite-size scaling behavior. The data points show a 
strong curvature in the log- log plot, so we first try an 
ansatz, including a single correction term, as 



C max (L) = C L a / v ll + ai L 



(16) 



with L being the stretched lattice length. Finite values 
for the specific heat exponent a do not lead to proper 
fits. Only for d — 7 the fit procedure converges to a/v = 
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FIG. 10: (color online) Collapsed susceptibility in d — 6 for 
several system sizes. The gray shaded area is the cumulative 
error of all data points, connected through c-splines, using 
v = 0.51(5), 7 = 1.07(5) and = 7.7859 
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FIG. 11: (color online) Specific heat in d = 5 for several 
system sizes. The lines are guides to the eyes only. 



-0.0016 ± 0.089, else it quits at a/v = 10~ 30 , leaving 
the other values unchanged. Peak heights and curves are 
shown in Fig. [lousing the values listed in Tab. lIIII Since 
the data of the specific heat at the peaks for the larger 
lattices are relatively rugged and the flanks are not, we 
had a look at the slopes of the specific heat. Interestingly, 
it seems that the minimum of the slope lies right at the 
critical point. But we did not elaborate on that point. 



d 


Co 


a/v 


Oi 


e 


5 


6.16(83) 


(fixed) 


-1.13(8) 


-0.34(9) 


6 


3.95(18) 


O(fixed) 


-1.17(6) 


-0.67(8) 


7 


3.1(9) 


-0.00(9) 


-1.6(4) 


-1.21(47) 



TABLE III: Results of the fits of the peak-height of the spe- 
cific heat according Eq. 



7 



3 - 



2.5 



o 



1.5 



- 


1 jy>^ 




/ HEH D=5 




/ I — Ah D=6 


II./ 1 


h-*-h D=7 
. . i . . . . i 



16 



FIG. 12: (color online) Specific heat peaks for d = 5, 6, 7 with 
fits according to Eq. (|16[) on log-log scale. 



To check our scaling assumptions we performed a data 
collapse. 

The obtained values lead to a good collapse for d = 5. 
For d = 6, 7 the data of the smaller system sizes do not 
collapse very well, but still OK. Using slightly changed 
parameters gives visually a better result, as shown in 
Fig. Q21 Nevertheless, the fit concerns only the correc- 
tions to scaling, and the improvement of the fit concerns 
the region away from the critical point. Therefore our 
conclusions are not changed. 
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FIG. 13: (color online) Collapse of the specific heat in d — 6 
(main plot) using e = —0.66 and in the inset d = 7 using 
e = -0.23. 



h c 
P 
7 

Q 
V 



d = 3 d = 4 d = 5 d = 6 d = 7 MF 

2.28(1) 4.18(1) 6.0157(10) 7.7859(1) 9.4889(1) 2d 1 /§ 

0.017(5) 0.13(5) 0.27(1) 0.50(1) 0.50(3) 1/2 

1.98(7) 1.57(10) 1.42(7) 1.07(3) 1.06(7) 1 



1.37(9) 0.78(10) 0.626(10) 0.51(5) 0.49(2) 1/2 



TABLE IV: Numerical results of the ground-state calculations 
for d = 34, 5, 6, 7 (data for d = 3 was taken from Refs^^, 
data for d = 4 was taken from Ref.— ) and the mean-field 
predictions (MF) for comparison. 



From those GSs we have calculated several thermody- 
namic quantities to determine the critical exponents that 
govern the FSS behaviour of the considered observables. 
A summary of the values thus obtained is given in Tab. 

m 

In this way, we have verified that the upper critical 
dimension of the RFIM is d u = 6 and that it exhibits 
mean-field critical behaviour at d u . In particular we ob- 
tained strong evidence for the specific heat of the RFIM 
to converge towards a constant for infinite systems sizes 
and found a = 0. Arguing from the extremely good data 
collapses in d — 6, a valid form of the logarithmic correc- 
tions of the disorder finite-size scaling for the short range 
RFIM has been found and quantified. 

For d = 7 we also find mean-field exponents when 
rescaling the system size according to Eq. ©, using a 
corrected correlation-length exponent for FSS analysis as 
^ (7) = ^mf d u /d. 

The results in d = 5 closed the gap to former simu- 
lations. Including them, a monotonic behaviour of the 
critical exponents as a function of the dimension can be 
seen. 

Finally, since the mean-field values for the exponents 
hold in d = 6, 7, the exponents do fulfill the Rushbrookc 
equation a + 2/3 + 7 = 2 exactly. For d = 5we also get 
a (5) + 2/3(5) + 7 (5) _ 2.0. 

To summarize, our results for high dimensions together 
with the results obtained earlier for d = 3,4 result in 
a rather complete picture of the RFIM order-disorder 
phase transition. 
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IV. CONCLUSIONS 



We have used a mapping to the max flow problem to 
calculate exact GSs of the RFIM at different disorder 
strenghts near the critical point in d = 5, 6, 7 dimensions. 
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